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We investigate the temperature dependence of optical absorption spectra of one-dimensional (ID) 
and two-dimensional (2D) Mott insulators by using an effective model in the strong-coupling limit 
of a half-filed Hubbard model. In the numerically exact diagonalization calculations on finite-size 
clusters, we find that in ID the energy position of the absorption edge is almost independent of 
temperature, while in 2D the edge position shifts to lower energy with increasing temperature. The 
different temperature dependence between ID and 2D is attributed to the difference of the coupling 
of the charge and spin degrees of freedom. The implications of the results on experiments are 
discussed in terms of the dimensionality dependence. 

PACS numbers: 78.30.Hv, 71.10.Fd, 78.20.Bh 



I. INTRODUCTION 

The charge gap in Mott insulators is a consequence of 
a strong electron correlation represented by a large on- 
site Coulomb interaction. The correlation induces novel 
phenomena in terms of the interplay of charge and spin 
degrees of freedom.^ In one-dimensional (ID) Mott insu- 
lators, two particles created by photoexcitation, i.e., an 
unoccupied site and a doubly occupied site of electrons, 
can move inside the system without being disturbed by 
surrounding spins in the background. This is a mani- 
festation of a separation of the charge and spin degrees 
of freedom, called the spin-charge separation. In two- 
dimensional (2D) Mott insulators, on the other hand, the 
two particles are expected not to be free from the spin 
degree of freedom, because the propagation of a carrier 
is known to induce a spin cloud around the carrier as 
a consequence of the misaligned spins along the carrier- 
hopping paths. Such an interplay of spin and charge 
is also one of the main subjects of the study of high- 
temperature superconductivity. 

The nature of the two particles in the photoexcited 
states of the Mott insulators is obtained by examining the 
linear susceptibility with respect to the applied electric 
field, which provides information on the dipole-allowed 
states with odd parity among the photoexcited states. In 
addition to the linear susceptibility, the third-order non- 
linear optical susceptibility is useful to detect not only 
the odd-parity states but also the dipole-forbidden states 
with even parity. Recently, large values of the nonlin- 
ear optical susceptibility have been reported for ID Mott 
insulators of copper oxides and halogen-bridged nickel 
compounds.^ ^ Analyses of the susceptibilities have sug- 
gested that odd- and even-parity states are nearly de- 
generate with a large transition dipole moment between 
them. Theoretically, the nonlinear susceptibility in ID 
Mott insulators has been examined by employing the nu- 
merically exact diagonalization technique for small clus- 
ters of the Hubbard model at half filling.^ It has been 
shown that odd- and even-parity states are almost degen- 
erate in the same energy region and that the degeneracy 



is due to the spin-charge separation and strong on-site 
Coulomb interaction.- 

In the 2D Mott insulators of copper oxide such as 
Sr2Cu02Cl2, the third-order nonlinear optical suscepti- 
bility has been reported to be one order of magnitude 
smaller than that in IDi^ Theoretical calculations based 
on the (extended) Hubbard model have shown such a 
dimensionality dependence!^ In addition to this, from 
the numerical investigation of the photoexcited states, it 
has been found that the edge of the distribution of the 
even-parity states is located lower in energy than that 
of the odd-parity states in 2D, in contrast to the almost 
degenerate behavior in IDpl The origin is attributed to 
the presence of the exchange interaction between spins, 
implying the important role of the spin degree of freedom 
in the photoexcited states. 

So far the dimensionality dependence of the photoex- 
cited states in the Mott insulators was theoretically ex- 
amined at zero temperature. With increasing tempera- 
ture, magnetic excitations are thermally created in the 
spin background irrespective of dimensionality. There- 
fore, it is not obvious whether the phtoexcited states are 
dependent on dimensionality even at finite temperatures. 
In this study, we numerically examine the temperature 
dependence of optical absorption spectra in a half-filled 
Hubbard model with large Coulomb interaction, and clar- 
ify its dimensionality dependence. It is found that there 
are qualitatively different behaviors between ID and 2D 
Mott insulators: in ID the energy position of the absorp- 
tion edge is almost independent of temperature as long 
as the temperature is below the exchange interaction be- 
tween spins, while in 2D the edge position shifts to lower 
energy with increasing temperature. Such a contrasting 
behavior is caused by the difference of the coupling of 
charge and spin that persists as long as short-range mag- 
netic interactions are present. The shift of the edge in 2D 
is consistent with experimental data in Sr2Cu02Cl2i^ 
The implications of our results to experiments will be 
discussed in terms of the dimensionality dependence and 
the effect of phonon scattering. 

The rest of this paper is organized as follows. We intro- 
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duce an effective Hamiltonian of the half-filled Hubbard 
model in the strong-cupling limit, and show outlines of 
the procedure for caluculating optical absorption spectra 
in Sec. II. In Sec. Ill, caluculated results of the tempera- 
ture dependence of optical absorption spectra in ID and 
2D sysytems are presents. The distributions of photoex- 
cited states are discussed in terms of the effect of the spin 
degree of freedom. The summary is given in Sec. IV. 



II. MODEL AND METHOD 

The Hubbard model is given by i?Hub ^ Ht + Hu 
with Ht = -tY.i,s.A(^\,<y(^i+S:<y + H.c.) and Hu = 
ni^^riij, where c\ ^ is the creation operator of an 
electron with spin cr at site i, ni_c = c\ ^Ci^a, i + S stands 
for the set of nearest neighbor (NN) sites of i, t is the 
NN hopping integral, and U is the on-site Coulomb in- 
teraction. 

In the strong coupling limit {U ^ t), there is no dou- 
bly occupied site in the initial states of optical absorp- 
tion at low temperatures (/cbT ^ U). (Hereafter, the 
Boltzmann factor is taken to be fee = !■) In this case, 
the initial states are described by the Heisenberg model 
given by Hq = Jj^i.si^i ' ^^+^ ^ fMnjA), where Si is 
the spin operator with S* = 1/2 at site z, rii = + '^i,!, 
and J = a£-IU. 

The photoexcited states, which are the final states of 
optical absorption, have both one doubly occupied site 
and one vacant site. An effective Hamiltonian describ- 
ing the photoexcited states is obtained by restricting the 
Hilbert spaces to a subspace with one doubly occupied 
site. By performing the second order perturbation with 
respect to the hopping term Ht in the Hubbard Hamil- 
tonian, the effective Hamiltonian is given by 
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where Hq, Hi, and H2 are projection operators onto the 
Hilbert space with zero, one, and two doubly occupied 
sites, respectively. In addition to them, we introduce 
an attractive interaction between a doubly occupied and 
vacant sites. This interaction is described by 
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The optical absorption spectrum at temperature T and 
photon energy uj is given by 



£2 (-,/?) = I (l-e-'^-)^e-^^^ 
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FIG. 1: The temperature dependence of the optical absorp- 
tion spectra from T = to T = 0.8t — 2.0J in a 16-site 
ID chain with U = Wt (J = OAt) and V = 1.5t. The delta- 
functions are broadened by a Lorentzian with a width of 0.05t. 



where p = 1/T, the initial state |/) is an eigenstate of Hq 
with energy Ej, and the final state \F) is an eigenstate of 
Hcs with energy Ep. The £/ represents spectral weight 
from each initial state |/). The energy Es is measured 
from the ground-state energy of Hq as Es — Ej — Eq, 
and Z = J2i e~^^^ . The is the x component of the 
current operator to second order of t. It is described hyi^ 

J- = if ^ + i_ J2 (s^-s'J 



The creation and annihilation operators, cj ^ and Ci^a, are 
projected onto the subspace with either zero or one dou- 
bly occupied site. The and S'^ are the x components 
of vectors connecting NN sites. 

We perform numerical calculations of the optical ab- 
sorption spectrum for 16-, 20-, and 24-site ID periodic 
chains, and \/20 x \/20- and \/26 x \/26- square clusters 
with periodic boundary conditions. For the 16-site chain, 
we calculate the optical absorption spectra by using all 
the initial states |/) and final states |-F) obtained by per- 
forming full diagonalizations of Hq and -ffcff + Hy- For 
the 20-site clusters, all the initial states |/) are obtained 
by diagonalizing Hq, but £/ for each |/) is calculated by 
the Lanczos technique. The Lanczos technique is also 
used for obtaining |/) in the case of the 24- and 26-site 
clusters. In this case, we obtain the low-energy eigen- 
states in the energy region of < Es < 4J as many as 
possible. 



III. RESULTS AND DISCUSSIONS 

Figure n shows £2 for the 16-site chain with U/t = 10 
{J/t = 0.4) and V/t = 1.5. At T = 0, there are only 
three peaks in the energy region oi uj/t < 10 because 
of finite-size system. With increasing temperature, the 
number of peaks increases according to the increase of 
the initial states contributing to the spectra. To be em- 
phasized here is that the position of absorption edges at 
around 6.5t is almost unchanged as long as the temper- 
ature is less than ^ 2J. This feature is independent of 
the size of the chains. The absence of the shift of the 
edge position is also seen in analytical calculationsii of 
the optical conductivity at T = and high temperatures 
(T ^ J). Note that at very high temperatures (T ^ t) 
the edge shifts to lower energy side accompanied by the 
broadening the spectrum. 
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FIG. 2: (a) ej in a 16-site ID chain with U = lOt (J = OAt) 
and V = 1.5t from various initial states |/) with energy Es 
and momentum k in units of tt. The ground-state energy is 
set to be Es = 0. (b) The size dependence of absorption- 
edge positions in a 16-site ID chains with J — OAt and J = 
O.OOlt {U — 400t). The circle, triangle, and square represent 
the positions from \I) with momentum A; = 0, 7r/2, and tt, 
respectively. The solid and open symbols denote the lowest- 
energy and next-lowest-energy peaks, respectively, for each 
1 7) as shown in (a). The solid and dotted lines are obtained 
by fitting them linearly to the symbols for J = OAt and J = 
0.00 If, respectively. 



FIG. 3: The temperature dependence of the optical absorp- 
tion spectra from T = to T = 0.8t — 2.0 J in (a) 26-site and 
(b) 20-site 2D square clusters with U = lOt (J = OAt) and 
V = 0. The delta-functions are broadened by a Lorentzian 
with a width of 0.05t. 



FIG. 4: £/ in (a) 26-site and (b) 20-site 2D square clusters 
with U = lot (J — OAt) and V — from various initial states 
|7) with energy Es and momentum k in units of n. The 
delta-functions are broadened by a Lorentzian with a width 
of 0.05*. 



In order to see the position of the edge more exphc- 
itly, we examine the contribution of each initial state |/) 
to the edge. Figuare[2Ia) shows decomposed spectra £/ 
from several low-energy initial state |/) that correspond 
to two-spinon excitations in the Heisenberg model. In 
FiglUa), J = OAt, and Es and k in each panel repre- 
sent the energy and the total momentum respectively of 
the corresponding initial state. The edge positions are 
found to be not much dependent on the initial states. 
To clarify the positions in the thermodynamic limit, we 
perform the size scaling of the positions of the lowest- 
energy and next-lowest-energy peaks for representative 
three total-momentum spaces {k = 0, 7r/2, and tt) in 
Fig. Elb). In — > cxD, N being system size, the lowest- 
energy peaks (solid symbols) converge to nearly the same 
position within a range of 0.03t irrespective of J/t (solid 
line: J/t = 0.4, and dotted line: J/t = 0.001). The 
next-lowest peaks also converge close to the lowest ones. 
We do not find any initial states that induce absorption 
edges lower in energy than that from the ground state. 
This is the reason why the absorption-edge position is 
independent of temperature as seen in Fig. ^ 

What is physics behind such little initial-state depen- 
dence of the edge position? At first sight, the positions 
seem to be dependent on the initial states because their 
energies are different from each other. However, the edge 
positions are actually almost independent of the initial 
states as seen in Fig. |21 (a). This means that the en- 
ergy variation among the final states which constitute 
the absorption edges is almost equivalent to that among 
the initial states. This is easily understood if we con- 
sider the fact that in the strong coupling limit the wave- 
functions in the ID Hubbard model can be given by the 
product of the charge part and the spin part, the latter 
of which is the wavefunction of the Heisenberg model. 
This leads to the same energy shift in the spin part for 
both the initial and final states. In other words, the dif- 
ference of the initial states as well as that of the final one 
comes from the different spin states but charge state is 
the same. Therefore, we can say that underlying physics 
behind no temperature dependence of the edge position 
is the separation of spin and charge degrees of freedom. 

Let us move to 2D systems. Figure O shows £2 for the 



%/26 X V26- and %/20 x ^/20-site clusters with U/t = 10 
{J/t = 0.4) and V/t = 0. At T = 0, spectral weights 
are predominantly concentrated on the absorption-edge 
region, although the spectral shapes slightly depend on 
the size of systems. The states near the edges may be de- 
scribed as excitonic states caused by antiferromagnetic 
interactions. With increasing temperature, spectral 
weights at absorption edge slightly sift to the lower en- 
ergy side. Other clusters such as 4 x 4 and vTS x vTS also 
show similar temperature dependence (not shown) . This 
dependence is in contrast with the ID systems discussed 
above. We decompose £2 into £/'s, which are shown in 
Fig. El for several initial states. The initial states shown 
are mainly the states constructing a spin- wave dispersion 
with 5 = 1 in the 2D Heisenberg model. In contrast to 
Fig. Ela), absorption-edge positions are found to be de- 
pendent on the initial states. Since in 2D the motion 
of photocarriers created in the final states induces mag- 
netic excitations unlike in ID, the energy variation among 
the final states are different from that among the initial 
states. This could be the origin of the shift of the edge 
positions. Therefore, it may be plausible to say that the 
coupling of spin and charge causes the temperature de- 
pendence of the edge positions. From Fig.01 we can also 
notice that sharp peaks at the edges (excitonic peaks) 
emerge depending on the initial states. By examining 
spin correlation functions of the initial states, we found 
that the excitonic peaks are pronounced in the absorp- 
tions from the initial states in which antiferromagnetic 
spin correlation remains to be long-range. This verifies 
an picture that the excitonic states are of magnetic origin 
and magnetically ordered states are necessary to realize 
them. 

Finally let us discuss the implications of the results to 
optical absorption experiments in Mott insulators. The 
temperature dependence of optical absorption has been 
reported for a 2D Mott insulator Sr2Cu02Cl2&^ The 
data show the shift of the edge toward lower-energy side 
with increasing temperature. This seems to be consistent 
with our results in 2D Hubbard model. However, the 
shift has been explained by taking into account phonon 
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scattering?^ Since we have not introduced phonon contri- 
butions in our model, we can not say the definite conclu- 
sion about the origin of the shift. From our results, we 
consider that not only the coupling with phonons but also 
with magnons^ may be necessary for precise description 
of the temperature dependence of the optical absorption. 
In ID, shifts of the edge positions similar to 2D systems 
have been reported. Since these data are different from 
our results, coupling with phonon should be examined. 

IV. SUMMARY 

In summary, we have examined the temperature de- 
pendence of optical absorption spectra of a half-filled ID 
and 2D Hubbard models in the strong-coupling limit. By 
employing the numerically exact diagonalization calcula- 
tions on finite-size clusters, we have found that in ID 
the energy position of the absorption edge is almost in- 



dependent of temperature, while in 2D the edge position 
shifts to lower energy with increasing temperature. The 
different temperature dependence between ID and 2D is 
attributed to the difference of the coupling of the charge 
and spin degrees of freedom. The analysis of existing 
experimental data is left until the effect of the electron- 
phonon coupling on absorption spectra is examined. 
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